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Using numerical holography, we study the collision of a planar sheet of energy with a bounded 
localized distribution of energy. The collision, which mimics proton-nucleus collisions, produces a 
localized lump of debris with transverse size R ~ 1/Teff with Tgff the effective temperature, and has 
large gradients and large transverse flow. Nevertheless, the postcollision evolution is well described 
by viscous hydrodynamics. Our results bolster the notion that debris produced in proton-nucleus 
collisions may be modeled using hydrodynamics. 


Introduction .—Data on recent light-heavy ion colli¬ 
sions, including proton-nucleus collisions, indicate the 
presence of collective flow nHi. While alternative mech¬ 
anisms for the origin of the flow have been proposed [5], 
the data are consistent with hydrodynamic evolution [i- 
[8]. However, at experimentally accessible energies micro¬ 
scopic scales are likely not too different than the system 
size. This raises several questions. Is it theoretically con¬ 
sistent to apply hydrodynamics to systems whose size is 
of order microscopic scales? Is it consistent to neglect 
nonhydrodynamic degrees of freedom when gradients are 
large? What is the size of the smallest drop of liquid? 

Since microscopic length and time scales typically de¬ 
crease in the limit of strong coupling, it is natural to ex¬ 
pect the domain of utility of hydrodynamics to be maxi¬ 
mized at strong coupling. However, strongly coupled dy¬ 
namics in QCD are notoriously difficult to study. Holo¬ 
graphic duality [9] maps the dynamics of certain strongly 
coupled non-Abelian gauge theories onto the dynamics of 
classical gravity in one higher dimension. The process of 
quark-gluon plasma formation is mapped onto gravita¬ 
tional collapse and black hole formation with the ring 
down of the black hole encoding the relaxation of the 
plasma to a hydrodynamic description. Since gravita¬ 
tional collapse can be studied numerically, holography 
provides a unique arena to study all stages of evolution 
— from far-from-equilibrium dynamics to hydrodynam¬ 
ics — in a microscopically complete and controlled set¬ 
ting. This has prompted much interest in using holo¬ 
graphic theories as toy models of real quark-gluon plasma 
(for a review see [lO]). The simplest theory with a dual 
holographic description is A? = 4 supersymmetric Yang- 
Mills theory (SYM), which is dual to gravity in AdSs. 

A simple model of quark-gluon plasma production is 
the collision of gravitational shock waves pni2Q] . which 
can result in the formation of a black hole. The profile 
of the gravitational waves encodes the profile of collid¬ 
ing shock waves in the dual field theory. In this Letter 
we report on shock wave collisions in SYM, which super¬ 
ficially at least, resemble proton-nucleus collisions, with 
the “nucleus” represented as a planar shock wave and 
the “proton” represented by shock wave localized in the 
transverse directions. We use quotes here to emphasize 


to the reader that the objects we are colliding are not 
bound states in QCD, but rather caricatures in SYM. 
The precollision geometry contains a trapped surface and 
the collision results in the formation of a black brane. We 
numerically solve the full 5D Einstein equations for the 
geometry after the collision and report on the evolution 
of the SYM stress tensor and test the validity of 
hydrodynamics. 

In strongly coupled SYM the hydrodynamic gradient 
expansion of is an expansion in powers of l/(ffeff) 
with £ the characteristic scale over which varies and 
Teff the effective temperature 1211122 ]. Hence in strongly 
coupled SYM the microscopic scale 1/Teff plays the role 
of the mean free path. We focus on the low energy limit, 
where Teff is small and 1/Teff is large, and demonstrate 
that collisions can result in the formation of droplets of 
liquid with sizes as small as i? ^ 1/Teff. Our results 
bolster the notion that debris produced in proton-nucleus 
collisions may be modeled using hydrodynamics. 

Gravitational formulation .— We construct initial data 
for Einstein’s equations by superimposing the metric of 
gravitational shock waves moving in the Hz directions at 
the speed of light. In Eefferman-Graham coordinates the 
metric of a single shock moving in the Hz direction is 


ds^ = [—dt^ H dx^ H ^ H h±{x_\_, Zz^,r) dz^] , (1) 

where x = {x, y, zj, x± = {x, y}, z^ = z ^t, and 


h±{x^ 





Shjk/r) 

^ 2^2 


( 2 ) 


The boundary of the spacetime lies at r = oo. The metric 
0 is an exact solution to Einstein’s equations for any 
choice of H± HD Eg. This geometry represents a state 
in the dual SYM theory with stress-energy tensor [29] 


^00 _ J.ZZ ^ ^J.0z _ H±(^X±,Z^), (3) 

(and all other components vanishing), where H± is the 
transverse Eourier transform of H±. We choose the waves 
moving in the Hz and — z directions to represent the “pro¬ 
ton” and “nucleus” respectively. 
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A simple choice of shock profiles is 

H±{xj_,z^) = (4) 

where is a smeared delta function with width w. 

On the boundary the longitudinally integrated energy 
density per unit area is Note that under a 

boost in the 2 ; direction (and in the w ^ 0 limit) only 
the normalization of change. We may therefore work 
in the frame in which max(/i+) = max(/i_). For nu¬ 
merical convenience we chose with 

<7 = 3, and = 1. Hence, the “proton” stress is 

localized about z = Xj_ = 0, and the “nucleus” stress 
is localized about z = —t and is translationally invari¬ 
ant in the transverse plane. Our choice of energy scale 
/i±{x±=0) = 1 fixes units in the results presented be¬ 
low. For the smeared longitudinal 6 function we use a 
Gaussian 6w(z) = . ^ with width w = 0.375. 

Note that evolution inside the future light cone of planar 
shock collisions with width re = 0.375 well approximates 
that of the 6 function limit [24] . 

For early times, t —re, the profiles H± have negli¬ 
gible overlap and the precollision geometry can be con¬ 
structed from § by replacing the last term with the 
sum of corresponding terms from left and right moving 
shocks. The resulting metric satisfies Einstein’s equa¬ 
tions, at early times, up to exponentially small errors. 

To evolve the precollision geometry forward in time 
we use the characteristic formulation of gravitational dy¬ 
namics in asymptotically AdS spacetimes discussed in 
detail in [25]. Our metric ansatz reads 

ds^ = r‘^ g^y{x^r) dx^dx^-\-2 dr dt ^ (5) 

with Greek indices denoting spacetime boundary coor¬ 
dinates, x^ = {t^x^y^z). Near the boundary, = 
/r^ + 0(l/r^). The subleading coefficients g^^J 
determine the SYM stress tensor, 

+ (6) 

To generate initial data for our characteristic evolu¬ 
tion, we numerically transform the precollision metric in 
Fefferman-Graham coordinates to the metric ansatz (§ 
m- We periodically compactify spatial directions with 
transverse size = Ly = 24 and longitudinal length 
Lz = 9. We begin time evolution at t = —1.5 and evolve 
to t = 2.5. Time evolution is performed using a spectral 
grid of size Nx = Ny = 39, Nz = 155 and Nr = 48. 

Results. — In Fig. [l] we plot the energy density 
(top) and momentum density (bottom) in the plane 
^ = 0 at several values of time. (Note that the plots are 
rotationally invariant about the 2 ; axis.) The color scaling 
in the lower plots denotes |T^*| and the flow lines indicate 
the direction of At time t = —1.125 the system con¬ 
sists of a planar sheet of energy (the “nucleus”) localized 


at 2 ; = 1.125 and a localized lump of energy (the “pro¬ 
ton”) centered at 2 ; = —1.125, x = 0, with transverse 
width cr = 3. These two distributions of energy move 
towards each other at the speed of light and collide at 
t = 2 : = 0. Far away from x = 0, where there is little over¬ 
lap between the “proton” and the “nucleus,” the collision 
has little effect on the future evolution of the “nucleus.” 
Indeed, at x = ±12 the “nucleus” simply continues to 
propagate at the speed of light with little change in the 
profile of the energy or momentum densities. However, 
at \x\ ^ a both the “proton” and “nucleus” are dramat¬ 
ically altered by the collision event; the collision results 
in the “proton” punching a hole of transverse size ~ a 
in the “nucleus.” The resulting produced debris lies in¬ 
side the light cone and subsequently expands both in the 
transverse and longitudinal directions. As we argue be¬ 
low, it is in the region \x\ ^ a — set by the transverse 
size of the “proton” — that the system begins to behave 
hydrodynamically after the collision. 

At time and length scales ^ than the microscopic scale 
1/Teff, nonhydro dynamic degrees of freedom must relax, 
the evolution of must be governed by hydrodynam¬ 
ics, and the dual black hole geometry must be governed 
by fluid/gravity duality [2TJ [22]. The salient hydrody¬ 
namic variables are the fluid velocity and the proper 
energy e, which near equilibrium are related to Teff by 

= (7) 

To obtain e and we extract the eigenvalues p(^x) and 
associated eigenvectors of 

e^A) = P(A) , (8) 

with no sum over A implied. Near equilibrium has one 
timelike eigenvector, and three spacelike eigenvec¬ 

tors, e^.y The temporal eigenvalue is simply the proper 
energy, e = —P(o), with the local fluid velocity the asso¬ 
ciated eigenvector, We choose normalization 

= — 1 with > 0. The spatial eigenvalues 
represent local pressures in the fluid rest frame. With 
the complete spacetime dependence of the flow field 
and energy density e determined from the exact stress, 
we then construct the hydrodynamic approximation to 
the stress tensor asing the constitutive relations of 

2 nd conformal hydrodynamics with transport coef¬ 

ficients determined by fluid/gravity duality [2TJ[22]. 

In Fig. [^ we plot the stress tensor components 
and , and their hydrodynamic approximations at cc = 
0 as a function of time. At this point u = 0(10“^), 
the stress tensor is approximately diagonal, and the local 
fluid pressures are well approximated by and 
The pressures increase dramatically during the collision, 
reflecting a system which is highly anisotropic and far 
from equilibrium. However, after time t = thydro ^ 1-2 
there is a qualitative change in the behavior of and 
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FIG. 1: The energy density (top) and momentum density |T°*| (bottom), at four different times, in the plane y = 0. 
Streamlines in the lower plots denote the direction of T°h The stress is rotationally invariant about the z axis. At time 
t = —1.125 the “nucleus” is at z = 1.125 and the “proton” is at z = —1.125. These two distributions of energy move towards 
each other at the speed of light and collide at t = z = 0. Note that the color scaling for the energy at t = 0 is off scale. Far 
away from x = 0, where there is little overlap between the “proton” and the “nucleus,” the collision has little effect on the 
future evolution of the “nucleus.” However, near x = 0 the “proton” punches a hole in the “nucleus.” The resulting debris has 
transverse size of order that of the “proton,” lies inside the light cone and expands in the transverse and longitudinal directions. 




t 


FIG. 2: Left: Stress tensor components and at cc = 0 
as a function of time together with their hydrodynamic ap¬ 
proximation. Around t = 0 the system is highly anisotropic 
and far from equilibrium. Nevertheless, at this point in space, 
the system begins to evolve hydrodynamically at t ^ 1.2. 
Right: the hydrodynamic residual A at cc = 0 as a function 
of time. After t — 1.2, A < 0.13. 


; thereafter the stress slowly varies in time and slowly 
isotropizes. Likewise, after thydro the pressures are well 
described the hydrodynamic constitutive relations. 

To elaborate on the domain of utility of hydrodynamics 
we define the dimensionless residual measure 


with ATA*- = TA*- - and p = I Ei PH) = f the av- 

erage pressure. In the local fluid rest frame AT^^ AT^^^ = 
AT'^^ ATij. Note that AT^^ AT^^, is in general nonmono¬ 
tonic in time due to nonhydro dynamic modes, such as 
quasinormal modes, oscillating in time while decaying. 
Indeed, and in Fig. agree with the hydrody¬ 
namic constitutive relations at t ^ 0.65, but subsequently 
disagree until t > 1.2. Taking the max in ameliorates 
the effect of the oscillations and aids in avoiding the false 
identification of hydrodynamic evolution if AT^^AT^^y is 
temporarily small. If A C 1 at t,cc, then is well 
described by the hydrodynamic constitutive relations at 
X at time t and all later times. Also included in Fig. is 
a plot of A at cc = 0 as a function of time. During the 
collision A ~ 1 and is not described by hydrodynam¬ 
ics. However, after t ^ 0.8, A rapidly decays and after 
t = 1.2, A < 0.13. In what follows we use A < 0.2 as a 
litmus test for whether the evolution of is consistent 
with the hydrodynamic constitutive relations. 

In the remaining plots Figs. we restrict our 

attention to time t = 1.5 and to the region in the x—z 
plane in where A < 0.2. We explain the coloring of the 
different plots below. The stress in the omitted regions 
has A > 0.2 and is not well approximated by the hydro- 
dynamic constitutive relations. All plots in Fig. [^include 
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FIG. 3: From left to right: the proper energy e, the anisotropy A, the first and second order gradient measures 7Z(i) and 7^(2), 
and the transverse and longitudinal fluid velocities Vx and Vz respectively. All plots are shown at time t = 1.5 and are restricted 
to the domain A < 0.2 and are rotationally invariant about the z axis. The outer (inner) contours in all plots are A = 0.2 (0.15), 
with A < 0.2(0.15) everywhere inside. In the displayed region, which has transverse radius R ~ 1/Teff, hydrodynamics is a 
good description of the evolution of the stress. 


the same A = const, contours. The outer (inner) con¬ 
tour is A = 0.2(0.15), with A < 0.2(0.15) everywhere 
inside. The transverse radius of these contours is ^ 4, 
which is approximately that of the initial “proton” seen 
in Fig.[^ We therefore conclude that the collision results 
in the formation of a droplet of liquid with transverse 
size of order that of the “proton.” The color scaling in 
Fig-i a) denotes the proper energy e. In the plotted do¬ 
main the average e is ~ 0.4 and, via Q, the average T^ff 
is ^ 0.25. We therefore obtain the dimensionless measure 
of the transverse size of the drop of liquid, 

RTeff 1 . ( 10 ) 

Evidently, hydrodynamics works even when the system 
size is of order the microscopic scale 1/Teff. Additionally, 
note thydroTeff ^ 0.3. Similar hydrodynamization times 
were observed in the 1+1 dimensional flows of [Hiiiiin] 
with thydroTeff ^ 0.3 appearing as a lower bound. 

When R ^ 1/Teff it is natural to expect gradients to 
be large. In order to quantify the size of gradients we 
define the anisotropy function 

A = max(p(i^)/inin(p(i)). ( 11 ) 

In ideal hydrodynamics, where the stress is isotropic in 
the local rest frame, A = 1. Therefore, deviations of A 
from 1 in regions where A <C 1 must be due to gradient 
corrections to the hydrodynamic constitutive relations. 
In Fig. [^b) we plot A. In the displayed domain A ^ 4. 
Evidently, gradients are large and ideal hydrodynamics 
is not a good approximation. 

However, the observed anisotropy is almost entirely 
due to first order viscous effects alone. To see this, let 
be the order contribution to the viscous stress 
and define the dimensionless gradient measures 

( 12 ) 


In the local rest frame +:) = 7^(„) 

measures the size of the order gradients correc¬ 
tions relative to ideal hydrodynamics. In Eigs. [^c)-[^d) 
we plot and 7^(2)- Inside the A = 0.15 contour 
7^(1) ^ 0.85 and 7^(2) ^0.1. Hence, first order gradient 
corrections are as large as the ideal stress with the second 
order correction an order of magnitude smaller. Similar 
observations were made in [m [22] within the context 
boost invariant flow. 

Discussion .—Given that 1/Teff is the salient micro¬ 
scopic scale in strongly coupled plasma — akin to a mean 
free path at weak coupling — it is remarkable that hy¬ 
drodynamics can describe the evolution systems as small 
as T ^ 1/Teff. Simply put, even when R ^ 1/Teff and 
gradients are large, nonhydro dynamic modes decay and 
the hydrodynamic gradient expansion is well behaved. 

It should be noted that upon fixing units by setting 
(7 equal to the proton radius, a ^ 1 fm, in the c.m. 
frame our “proton” has energy ^ 20 GeV and the ef¬ 
fective temperature of the produced plasma is Teff ^ 200 
MeV. Given that proton-nucleus collisions at RHIC and 
the LHC have energies and temperatures greater than 
this, it is natural to expect TTeff to be larger in RHIC 
and LHC collisions than the simulated collision presented 
here. Given that hydrodynamics already works well when 
TTeff ^ 1, our observed extreme applicability of hydro¬ 
dynamics, both in terms of the system size and the size 
of gradients, bolsters the notion the debris produced in 
proton-nucleus collisions can be modeled with hydrody¬ 
namics. 

However, hydrodynamic simulations of tiny short lived 
systems with large gradients will likely be sensitive to the 
initial viscous stress, which Eigs.[^b)-|^c) demonstrate is 
large, and to the initial transverse fluid velocity. Indeed, 
as shown in Eig. ie)i:f), the transverse and longitu¬ 
dinal components of the fluid 3-velocity v = u/u^ are 
similar in magnitude. Simply put, large gradients result 
in large viscous stress and drive the rapid development 
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of transverse flow. These effects are likely key ingredi¬ 
ents in initial hydrodynamic data required for precision 
modeling of proton-nucleus collisions. 

It would be interesting to push our analysis further: 
how big are the smallest drops of liquid? Clearly i?Teff 
cannot be made arbitrarily small since the hydrody¬ 
namic gradient expansion is badly behaved when i?Teff ^ 
1. Additionally, in the gravitational description there 
should exist a critical energy below which no black 
hole is formed, and which for ^ = £’c + 0^ critical gravi¬ 
tational collapse occurs m- The absence of a black hole 
when E < Ec means that in the dual field theory the col- 
lisional debris will not evolve hydro dynamically at any 
future time. It would be interesting to study dynamics 
near Ec and the signatures of hydrodynamics turning off. 
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